Taxonomic distribution of metabolic functions in bacteria associated with Trichodesmium consortia.

IMPORTANCE: Colonies of the cyanobacteria Trichodesmium act as a biological hotspot for the usage and recycling of key resources such as C, N, P, and Fe within an otherwise oligotrophic environment. While Trichodesmium colonies are known to interact and support a unique community of algae and particle-associated microbes, our understanding of the taxa that populate these colonies and the gene functions they encode is still limited. Characterizing the taxa and adaptive strategies that influence consortium physiology and its concomitant biogeochemistry is critical in a future ocean predicted to have increasingly resource-depleted regions.

Microbial community transcriptional patterns vary in response to mesoscale forcing in the North Pacific Subtropical Gyre.

The physical and biological dynamics that influence phytoplankton communities in the oligotrophic ocean are complex, changing across broad temporal and spatial scales. Eukaryotic phytoplankton (e.g., diatoms), despite their relatively low abundance in oligotrophic waters, are responsible for a large component of the organic matter flux to the ocean interior. Mesoscale eddies can impact both microbial community structure and function, enhancing primary production and carbon export, but the mechanisms that underpin these dynamics are still poorly understood. Here, mesoscale eddy influences on the taxonomic diversity and expressed functional profiles of surface communities of microeukaryotes and particle-associated heterotrophic bacteria from the North Pacific Subtropical Gyre were assessed over 2 years (spring 2016 and summer 2017). The taxonomic diversity of the microeukaryotes significantly differed by eddy polarity (cyclonic versus anticyclonic) and between sampling seasons/years and was significantly correlated with the taxonomic diversity of particle-associated heterotrophic bacteria. The expressed functional profile of these taxonomically distinct microeukaryotes varied consistently as a function of eddy polarity, with cyclones having a different expression pattern than anticyclones, and between sampling seasons/years. These data suggest that mesoscale forcing, and associated changes in biogeochemistry, could drive specific physiological responses in the resident microeukaryote community, independent of species composition.

Transcriptional response of Emiliania huxleyi under changing nutrient environments in the North Pacific Subtropical Gyre.

The widespread coccolithophore Emiliania huxleyi is an abundant oceanic phytoplankton, impacting the global cycling of carbon through both photosynthesis and calcification. Here, we examined the transcriptional responses of populations of E. huxleyi in the North Pacific Subtropical Gyre to shifts in the nutrient environment. Using a metatranscriptomic approach, nutrient-amended microcosm studies were used to track the global metabolism of E. huxleyi. The addition of nitrate led to significant changes in transcript abundance for gene pathways involved in nitrogen and phosphorus metabolism, with a decrease in the abundance of genes involved in the acquisition of nitrogen (e.g. N-transporters) and an increase in the abundance of genes associated with phosphate acquisition (e.g. phosphatases). Simultaneously, after the addition of nitrate, genes associated with calcification and genes unique to the diploid life stages of E. huxleyi significantly increased. These results suggest that nitrogen is a major driver of the physiological ecology of E. huxleyi in this system and further suggest that the addition of nitrate drives shifts in the dominant life-stage of the population. Together, these results underscore the importance of phenotypic plasticity to the success of E. huxleyi, a characteristic that likely underpins its ability to thrive across a variety of marine environments.

Transcriptional patterns of Emiliania huxleyi in the North Pacific Subtropical Gyre reveal the daily rhythms of its metabolic potential.

Emiliania huxleyi is a calcifying haptophyte, contributing to both the organic and inorganic marine carbon cycles. In marine ecosystems, light is a major driver of phytoplankton physiology and ultimately carbon flow through the ecosystem. Here, we analysed a Lagrangian time-series of metatranscriptomes collected in the North Pacific Subtropical Gyre (NPSG) to examine how in situ populations of E. huxleyi modulate gene expression over day-night transitions. Many E. huxleyi contigs had a diel expression pattern, with 61% of contigs clustering into modules with statistically significant diel periodicity. Contigs involved in processes that build up energy stores, like carbon fixation and lipid synthesis, peaked around dawn. In contrast, contigs involved in processes that released energy stores, like respiration and lipid degradation, peaked mid-day and towards dusk. These patterns suggest an orchestrated cycle of building, then consuming energy stores in E. huxleyi populations in the NPSG. Selected contigs related to the cell cycle also exhibited significant diel periodicity consistent with phased modulations of division observed in culture. Overall, these patterns of gene expression suggest a daily metabolic cascade that could contribute to both organic and inorganic carbon flow in this nutrient depleted ecosystem.

Transcriptional and Proteomic Choreography Under Phosphorus Deficiency and Re-supply in the N2 Fixing Cyanobacterium Trichodesmium erythraeum.

The N2 fixing cyanobacterium Trichodesmium is a critically important organism in oligotrophic marine ecosystems, supplying "new" nitrogen (N) to the otherwise N-poor tropical and subtropical regions where it occurs. Low concentrations of phosphorus (P) in these regions can constrain Trichodesmium distribution and N2 fixation rates. Physiological characterization of a single species in a mixed community can be challenging, and 'omic approaches are increasingly important tools for tracking nutritional physiology in a taxon-specific manner. As such, studies examining the dynamics of gene and protein markers of physiology (e.g., nutrient stress) are critical for the application and interpretation of such 'omic data in situ. Here we leveraged combined transcriptomics, proteomics, and enzyme activity assays to track the physiological response of Trichodesmium erythraeum IMS101 to P deficiency and subsequent P re-supply over 72 h of sampling. P deficiency resulted in differential gene expression, protein abundance, and enzyme activity that highlighted a synchronous shift in P physiology with increases in the transcripts and corresponding proteins for hydrolyzing organic phosphorus, taking up phosphate with higher affinity, and modulating intracellular P demand. After P deficiency was alleviated, gene expression of these biomarkers was reduced to replete levels within 4 h of P amendment. A number of these gene biomarkers were adjacent to putative pho boxes and their expression patterns were similar to a sphR response regulator. Protein products of the P deficiency biomarkers were slow to decline, with 84% of the original P deficient protein set still significantly differentially expressed after 72 h. Alkaline phosphatase activity tracked with proteins for this enzyme. With the rapid turnover time of transcripts, they appear to be good biomarkers of a P stress phenotype, whereas proteins, with a slower turnover time, may better reflect cellular activities. These results highlight the importance of validating and pairing transcriptome and proteome data that can be applied to physiological studies of key species in situ.

Transcriptional Shifts Highlight the Role of Nutrients in Harmful Brown Tide Dynamics.

Harmful algal blooms (HABs) threaten ecosystems and human health worldwide. Controlling nitrogen inputs to coastal waters is a common HAB management strategy, as nutrient concentrations often suggest coastal blooms are nitrogen-limited. However, defining best nutrient management practices is a long-standing challenge: in part, because of difficulties in directly tracking the nutritional physiology of harmful species in mixed communities. Using metatranscriptome sequencing and incubation experiments, we addressed this challenge by assaying the in situ physiological ecology of the ecosystem destructive alga, Aureococcus anophagefferens. Here we show that gene markers of phosphorus deficiency were expressed in situ, and modulated by the enrichment of phosphorus, which was consistent with the observed growth rate responses. These data demonstrate the importance of phosphorus in controlling brown-tide dynamics, suggesting that phosphorus, in addition to nitrogen, should be evaluated in the management and mitigation of these blooms. Given that nutrient concentrations alone were suggestive of a nitrogen-limited ecosystem, this study underscores the value of directly assaying harmful algae in situ for the development of management strategies.

Non-linear Physiology and Gene Expression Responses of Harmful Alga Heterosigma akashiwo to Rising CO2.

Heterosigma akashiwo is a raphidophyte known for forming ichthyotoxic blooms. In order to predict the potential impacts of rising CO2 on H. akashiwo it is necessary to understand the factors influencing growth rates over a range of CO2 concentrations. Here we examined the physiology and gene expression response of H. akashiwo to concentrations from 200 to 1000ppm CO2. Growth rate data were combined from this and previous studies and fit with a CO2 limitation-inhibition model that revealed an apparent growth optimum around 600-800ppm CO2. Physiological changes included a significant increase in C:N ratio at ∼800ppm CO2 and a significant decrease in hydrogen peroxide concentration at ∼1000ppm. Whole transcriptome sequencing of H. akashiwo revealed sharp distinctions in metabolic pathway gene expression between ∼600 and ∼800ppm CO2. Hierarchical clustering by co-expression identified groups of genes with significant correlations to CO2 and growth rate. Genes with significant differential expression with CO2 included carbon concentrating mechanism genes such as beta-carbonic anhydrases and a bicarbonate transporter, which may underpin shifts in physiology. Genes involved in cell motility were significantly changed by both elevated CO2 and growth rate, suggesting that future ocean conditions could modify swimming behavior in this species.

Periodic and coordinated gene expression between a diazotroph and its diatom host.

In the surface ocean, light fuels photosynthetic carbon fixation of phytoplankton, playing a critical role in ecosystem processes including carbon export to the deep sea. In oligotrophic oceans, diatom-diazotroph associations (DDAs) play a keystone role in ecosystem function because diazotrophs can provide otherwise scarce biologically available nitrogen to the diatom host, fueling growth and subsequent carbon sequestration. Despite their importance, relatively little is known about the nature of these associations in situ. Here we used metatranscriptomic sequencing of surface samples from the North Pacific Subtropical Gyre (NPSG) to reconstruct patterns of gene expression for the diazotrophic symbiont Richelia and we examined how these patterns were integrated with those of the diatom host over day-night transitions. Richelia exhibited significant diel signals for genes related to photosynthesis, N2 fixation, and resource acquisition, among other processes. N2 fixation genes were significantly co-expressed with host nitrogen uptake and metabolism, as well as potential genes involved in carbon transport, which may underpin the exchange of nitrogen and carbon within this association. Patterns of expression suggested cell division was integrated between the host and symbiont across the diel cycle. Collectively these data suggest that symbiont-host physiological ecology is strongly interconnected in the NPSG.

Transcriptional patterns identify resource controls on the diazotroph Trichodesmium in the Atlantic and Pacific oceans.

The N2-fixing cyanobacterium Trichodesmium is intensely studied because of the control this organism exerts over the cycling of carbon and nitrogen in the low nutrient ocean gyres. Although iron (Fe) and phosphorus (P) bioavailability are thought to be major drivers of Trichodesmium distributions and activities, identifying resource controls on Trichodesmium is challenging, as Fe and P are often organically complexed and their bioavailability to a single species in a mixed community is difficult to constrain. Further, Fe and P geochemistries are linked through the activities of metalloenzymes, such as the alkaline phosphatases (APs) PhoX and PhoA, which are used by microbes to access dissolved organic P (DOP). Here we identified significant correlations between Trichodesmium-specific transcriptional patterns in the North Atlantic (NASG) and North Pacific Subtropical Gyres (NPSG) and patterns in Fe and P biogeochemistry, with the relative enrichment of Fe stress markers in the NPSG, and P stress markers in the NASG. We also observed the differential enrichment of Fe-requiring PhoX transcripts in the NASG and Fe-insensitive PhoA transcripts in the NPSG, suggesting that metalloenzyme switching may be used to mitigate Fe limitation of DOP metabolism in Trichodesmium. This trait may underpin Trichodesmium success across disparate ecosystems.

Coordinated gene expression between Trichodesmium and its microbiome over day-night cycles in the North Pacific Subtropical Gyre.

Trichodesmium is a widespread, N2 fixing marine cyanobacterium that drives inputs of newly fixed nitrogen and carbon into the oligotrophic ecosystems where it occurs. Colonies of Trichodesmium ubiquitously occur with heterotrophic bacteria that make up a diverse microbiome, and interactions within this Trichodesmium holobiont could influence the fate of fixed carbon and nitrogen. Metatranscriptome sequencing was performed on Trichodesmium colonies collected during high-frequency Lagrangian sampling in the North Pacific Subtropical Gyre (NPSG) to identify possible interactions between the Trichodesmium host and microbiome over day-night cycles. Here we show significantly coordinated patterns of gene expression between host and microbiome, many of which had significant day-night periodicity. The functions of the co-expressed genes suggested a suite of interactions within the holobiont linked to key resources including nitrogen, carbon, and iron. Evidence of microbiome reliance on Trichodesmium-derived vitamin B12 was also detected in co-expression patterns, highlighting a dependency that could shape holobiont community structure. Collectively, these patterns of expression suggest that biotic interactions could influence colony cycling of resources like nitrogen and vitamin B12, and decouple activities, like N2 fixation, from typical abiotic drivers of Trichodesmium physiological ecology.

The impact of elevated CO2 on Prochlorococcus and microbial interactions with 'helper' bacterium Alteromonas.

Prochlorococcus is a globally important marine cyanobacterium that lacks the gene catalase and relies on 'helper' bacteria such as Alteromonas to remove reactive oxygen species. Increasing atmospheric CO2 decreases the need for carbon concentrating mechanisms and photorespiration in phytoplankton, potentially altering their metabolism and microbial interactions even when carbon is not limiting growth. Here, Prochlorococcus (VOL4, MIT9312) was co-cultured with Alteromonas (strain EZ55) under ambient (400 p.p.m.) and elevated CO2 (800 p.p.m.). Under elevated CO2, Prochlorococcus had a significantly longer lag phase and greater apparent die-offs after transfers suggesting an increase in oxidative stress. Whole-transcriptome analysis of Prochlorococcus revealed decreased expression of the carbon fixation operon, including carboxysome subunits, corresponding with significantly fewer carboxysome structures observed by electron microscopy. Prochlorococcus co-culture responsive gene 1 had significantly increased expression in elevated CO2, potentially indicating a shift in the microbial interaction. Transcriptome analysis of Alteromonas in co-culture with Prochlorococcus revealed decreased expression of the catalase gene, known to be critical in relieving oxidative stress in Prochlorococcus by removing hydrogen peroxide. The decrease in catalase gene expression was corroborated by a significant ~6-fold decrease in removal rates of hydrogen peroxide from co-cultures. These data suggest Prochlorococcus may be more vulnerable to oxidative stress under elevated CO2 in part from a decrease in ecosystem services provided by heterotrophs like Alteromonas. This work highlights the importance of considering microbial interactions in the context of a changing ocean.The ISME Journal advance online publication, 31 October 2017; doi:10.1038/ismej.2017.189.

Transcriptional response of the harmful raphidophyte Heterosigma akashiwo to nitrate and phosphate stress.

The marine eukaryotic alga Heterosigma akashiwo (Raphidophyceae) is known for forming ichthyotoxic harmful algal blooms (HABs). In the past 50 years, H. akashiwo blooms have increased, occurring globally in highly eutrophic coastal and estuarine systems. These systems often incur dramatic physicochemical changes, including macronutrient (nitrogen and phosphorus) enrichment and depletion, on short timescales. Here, H. akashiwo cultures grown under nutrient replete, low N and low P growth conditions were examined for changes in biochemical and physiological characteristics in concert with transcriptome sequencing to provide a mechanistic perspective on the metabolic processes involved in responding to N and P stress. There was a marked difference in the overall transcriptional pattern between low N and low P transcriptomes. Both nutrient stresses led to significant changes in the abundance of thousands of contigs related to a wide diversity of metabolic pathways, with limited overlap between the transcriptomic responses to low N and low P. Enriched contigs under low N included many related to nitrogen metabolism, acquisition, and transport. In addition, metabolic modules like photosynthesis and carbohydrate metabolism changed significantly under low N, coincident with treatment-specific changes in photosynthetic efficiency and particulate carbohydrate content. P-specific contigs responsible for P transport and organic P use were more enriched in the low P treatment than in the replete control and low N treatment. These results provide new insight into the genetic mechanisms that distinguish how this HAB species responds to these two common nutrient stresses, and the results can inform future field studies, linking transcriptional patterns to the physiological ecology of H. akashiwo in situ.

Conserved Transcriptional Responses to Nutrient Stress in Bloom-Forming Algae.

The concentration and composition of bioavailable nitrogen (N) and phosphorus (P) in the upper ocean shape eukaryotic phytoplankton communities and influence their physiological responses. Phytoplankton are known to exhibit similar physiological responses to limiting N and P conditions such as decreased growth rates, chlorosis, and increased assimilation of N and P. Are these responses similar at the molecular level across multiple species? To interrogate this question, five species from biogeochemically important, bloom-forming taxa (Bacillariophyta, Dinophyta, and Haptophyta) were grown under similar low N, low P, and replete nutrient conditions to identify transcriptional patterns and associated changes in biochemical pools related to N and P stress. Metabolic profiles, revealed through the transcriptomes of these taxa, clustered together based on species rather than nutrient stressor, suggesting that the global metabolic response to nutrient stresses was largely, but not exclusively, species-specific. Nutrient stress led to few transcriptional changes in the two dinoflagellates, consistent with other research. An orthologous group analysis examined functionally conserved (i.e., similarly changed) responses to nutrient stress and therefore focused on the diatom and haptophytes. Most conserved ortholog changes were specific to a single nutrient treatment, but a small number of orthologs were similarly changed under both N and P stress in 2 or more species. Many of these orthologs were related to photosynthesis and may represent generalized stress responses. A greater number of orthologs were conserved across more than one species under low P compared to low N. Screening the conserved orthologs for functions related to N and P metabolism revealed increased relative abundance of orthologs for nitrate, nitrite, ammonium, and amino acid transporters under N stress, and increased relative abundance of orthologs related to acquisition of inorganic and organic P substrates under P stress. Although the global transcriptional responses were dominated by species-specific changes, the analysis of conserved responses revealed functional similarities in resource acquisition pathways among different phytoplankton taxa. This overlap in nutrient stress responses observed among species may be useful for tracking the physiological ecology of phytoplankton field populations.

Diverse CO2-Induced Responses in Physiology and Gene Expression among Eukaryotic Phytoplankton.

With rising atmospheric CO2, phytoplankton face shifts in ocean chemistry including increased dissolved CO2 and acidification that will likely influence the relative competitive fitness of different phytoplankton taxa. Here we compared the physiological and gene expression responses of six species of phytoplankton including a diatom, a raphidophyte, two haptophytes, and two dinoflagellates to ambient (~400 ppm) and elevated (~800 ppm) CO2. Dinoflagellates had significantly slower growth rates and higher, yet variable, chlorophyll a per cell under elevated CO2. The other phytoplankton tended to have increased growth rates and/or decreased chlorophyll a per cell. Carbon and nitrogen partitioning of cells shifted under elevated CO2 in some species, indicating potential changes in energy fluxes due to changes in carbon concentrating mechanisms (CCM) or photorespiration. Consistent with these phenotypic changes, gene set enrichment analyses revealed shifts in energy, carbon and nitrogen metabolic pathways, though with limited overlap between species in the genes and pathways involved. Similarly, gene expression responses across species revealed few conserved CO2-responsive genes within CCM and photorespiration categories, and a survey of available transcriptomes found high diversity in biophysical CCM and photorespiration expressed gene complements between and within the four phyla represented by these species. The few genes that displayed similar responses to CO2 across phyla were from understudied gene families, making them targets for further research to uncover the mechanisms of phytoplankton acclimation to elevated CO2. These results underscore that eukaryotic phytoplankton have diverse gene complements and gene expression responses to CO2 perturbations and highlight the value of cross-phyla comparisons for identifying gene families that respond to environmental change.

Variable depth distribution of Trichodesmium clades in the North Pacific Ocean.

Populations of nitrogen-fixing cyanobacteria in the genus Trichodesmium are critical to ocean ecosystems, yet predicting patterns of Trichodesmium distribution and their role in ocean biogeochemistry is an ongoing challenge. This may, in part, be due to differences in the physiological ecology of Trichodesmium species, which are not typically considered independently in field studies. In this study, the abundance of the two dominant Trichodesmium clades (Clade I and Clade III) was investigated during a survey at Station ALOHA in the North Pacific Subtropical Gyre (NPSG) using a clade-specific qPCR approach. While Clade I dominated the Trichodesmium community, Clade III abundance was >50% in some NPSG samples, in contrast to the western North Atlantic where Clade III abundance was always <10%. Clade I populations were distributed down to depths >80 m, while Clade III populations were only observed in the mixed layer and found to be significantly correlated with depth and temperature. These data suggest active niche partitioning of Trichodesmium species from different clades, as has been observed in other cyanobacteria. Tracking the distribution and physiology of Trichodesmium spp. would contribute to better predictions of the physiological ecology of this biogeochemically important genus in the present and future ocean.

Microbial diversity within the Trichodesmium holobiont.

Nitrogen-fixing cyanobacteria in the genus Trichodesmium play a critical role in the productivity of the tropical and subtropical oligotrophic oceans. The ecological success of these populations is likely associated with the diverse microbial interactions occurring within the Trichodesmium holobiont, especially between Trichodesmium and heterotrophic bacterial epibionts. Yet, the composition of the Trichodesmium holobiont and the processes governing microbial assemblage are not well documented. Here, we used high-resolution 16S rDNA amplicon sequencing to examine the diversity of Trichodesmium and associated epibionts across different ocean regions and colony morphologies (puffs and rafts). Trichodesmium Clade I (i.e., T. thiebautii-like) dominated the colonies in all ocean basins regardless of morphology, although the Trichodesmium community structure significantly varied between morphologies in some regions. On average, Alphaproteobacteria (i.e., Thalassobius), Gammaproteobacteria (i.e., Pseudoalteromonas), Sphingobacteria (i.e., Microscilla and Vibrio) and Flavobacteria dominated the epibiont communities, but community composition and structure significantly differed between regions. Epibionts from the two colony morphologies were taxonomically and functionally distinct in the North Atlantic and North Pacific. These findings suggest that the colony types might define two distinct niches and that epibiont assemblage might be driven in part by selective processes, where epibionts are selected according to their influence on colony metabolism.

Functional group-specific traits drive phytoplankton dynamics in the oligotrophic ocean.

A diverse microbial assemblage in the ocean is responsible for nearly half of global primary production. It has been hypothesized and experimentally demonstrated that nutrient loading can stimulate blooms of large eukaryotic phytoplankton in oligotrophic systems. Although central to balancing biogeochemical models, knowledge of the metabolic traits that govern the dynamics of these bloom-forming phytoplankton is limited. We used eukaryotic metatranscriptomic techniques to identify the metabolic basis of functional group-specific traits that may drive the shift between net heterotrophy and autotrophy in the oligotrophic ocean. Replicated blooms were simulated by deep seawater (DSW) addition to mimic nutrient loading in the North Pacific Subtropical Gyre, and the transcriptional responses of phytoplankton functional groups were assayed. Responses of the diatom, haptophyte, and dinoflagellate functional groups in simulated blooms were unique, with diatoms and haptophytes significantly (95% confidence) shifting their quantitative metabolic fingerprint from the in situ condition, whereas dinoflagellates showed little response. Significantly differentially abundant genes identified the importance of colimitation by nutrients, metals, and vitamins in eukaryotic phytoplankton metabolism and bloom formation in this system. The variable transcript allocation ratio, used to quantify transcript reallocation following DSW amendment, differed for diatoms and haptophytes, reflecting the long-standing paradigm of phytoplankton r- and K-type growth strategies. Although the underlying metabolic potential of the large eukaryotic phytoplankton was consistently present, the lack of a bloom during the study period suggests a crucial dependence on physical and biogeochemical forcing, which are susceptible to alteration with changing climate.

Arsenate Resistance in the Unicellular Marine Diazotroph Crocosphaera watsonii.

The toxic arsenate ion can behave as a phosphate analog, and this can result in arsenate toxicity especially in areas with elevated arsenate to phosphate ratios like the surface waters of the ocean gyres. In these systems, cellular arsenate resistance strategies would allow phytoplankton to ameliorate the effects of arsenate transport into the cell. Despite the potential coupling between arsenate and phosphate cycling in oligotrophic marine waters, relatively little is known about arsenate resistance in the nitrogen-fixing marine cyanobacteria that are key components of the microbial community in low nutrient systems. The unicellular diazotroph, Crocosphaera watsonii WH8501, was able to grow at reduced rates with arsenate additions up to 30 nM, and estimated arsenate to phosphate ratios of 6:1. The genome of strain WH8501 contains homologs for arsA, arsH, arsB, and arsC, allowing for the reduction of arsenate to arsenite and the pumping of arsenite out of the cell. The short-term addition of arsenate to the growth medium had no effect on nitrogen fixation. However, arsenate addition did result in the up-regulation of the arsB gene with increasing arsenate concentrations, indicating the induction of the arsenate detoxification response. The arsB gene was also up-regulated by phosphorus stress in concert with a gene encoding the high-affinity phosphate binding protein pstS. Both genes were down-regulated when phosphate was re-fed to phosphorus-stressed cells. A field survey of surface water from the low phosphate western North Atlantic detected expression of C. watsoniiarsB, suggestive of the potential importance of arsenate resistance strategies in this and perhaps other systems.

Nutrient-regulated transcriptional responses in the brown tide-forming alga Aureococcus anophagefferens.

Long-SAGE (serial analysis of gene expression) was used to profile the transcriptome of the brown tide-forming alga, Aureococcus anophagefferens, under nutrient replete (control), and nitrogen (N) and phosphorus (P) deficiency to understand how this organism responds at the transcriptional level to varying nutrient conditions. This approach has aided A. anophagefferens genome annotation efforts and identified a suite of genes upregulated by N and P deficiency, some of which have known roles in nutrient metabolism. Genes upregulated under N deficiency include an ammonium transporter, an acetamidase/formamidase and two peptidases. This suggests an ability to utilize reduced N compounds and dissolved organic nitrogen, supporting the hypothesized importance of these N sources in A. anophagefferens bloom formation. There are also a broad suite of P-regulated genes, including an alkaline phosphatase, and two 5'-nucleotidases, suggesting A. anophagefferens may use dissolved organic phosphorus under low phosphate conditions. These N- and P-regulated genes may be important targets for exploring nutrient controls on bloom formation in field populations.

Parallel analyses of Alexandrium catenella cell concentrations and shellfish toxicity in the Puget Sound.

Alexandrium catenella is widespread in western North America and produces a suite of potent neurotoxins that cause paralytic shellfish poisoning (PSP) in humans and have deleterious impacts on public health and economic resources. There are seasonal PSP-related closures of recreational and commercial shellfisheries in the Puget Sound, but the factors that influence cell distribution, abundance, and relationship to paralytic shellfish toxins (PSTs) in this system are poorly described. Here, a quantitative PCR assay was used to detect A. catenella cells in parallel with state shellfish toxicity testing during the 2006 bloom season at 41 sites from April through October. Over 500,000 A. catenella cells liter(-1) were detected at several stations, with two main pulses of cells driving cell distribution, one in June and the other in August. PSTs over the closure limit of 80 mug of PST 100 per g of shellfish tissue were detected at 26 of the 41 sites. Comparison of cell numbers and PST data shows that shellfish toxicity is preceded by an increase in A. catenella cells in 71% of cases. However, cells were also observed in the absence of PSTs in shellfish, highlighting the complex relationship between A. catenella and the resulting shellfish toxicity. These data provide important information on the dynamics of A. catenella cells in the Puget Sound and are a first step toward assessing the utility of plankton monitoring to augment shellfish toxicity testing in this system.